Multi-octave supercontinuum from mid-IR filamentation in bulk 
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In supercontinuum generation, various propa- 
gation effects combine to produce a dramatic spec- 
tral broadening [ ] of intense ultrashort optical 
pulses with far reaching possibilities [ ]. Different 
applications place highly divergent and challenging 
demands on source characteristics such as spectral 
coverage from the ultraviolet (UV) 4] across 
the visible (VIS) to the near-infrared (NIR), and 
into the mid-infrared (MIR) [5, 6, 7, 8]. Shot-to- 
shot repeatability, high spectral energy density, an 
absence of complicated or non-deterministic pulse 
splitting are also essential for many applications. 
Here we present an "all in one" solution with the 
first supercontinuum in bulk covering the broad- 
est bandwidth from just above UV far into the 
MIR. The spectrum spans more than three oc- 
taves, carries high spectral energy density (3 pJ 
up to 10 nJ per nanometer), has high shot-to-shot 
reproducibility, and is carrier-to-envelope phase 
(CEP) stable. Our method, based on filamenta- 
tion [9] of a femtosecond MIR pulse in the anoma- 
lous dispersion regime, allows for a new class of 
simple and compact supercontinuum sources. 

Coherent ultra-broadband light sources are in demand 
for a large variety of applications. Modern imaging — 
especially optical coherence tomography [ ] — exploits 
broadband light sources to map biological materials with 
high spatial resolution. Multi-photon spectroscopy and 
microscopy [ ] demand focusable and coherent light at 
a variety of wavelengths to visualize cells or for chem- 
ical imaging of cellular components. Modern nonlinear 
spectroscopy [ ] benefits from time-resolved multi-color 
pump-probe capabilities for tracking the motion of va- 
lence electrons and nuclei on their natural femtosecond 
timescale. Breath analysis is a promising technique for 
early cancer detection demanding spectroscopic detection 
of volatiles which is achieved with high sensitivity in the 
molecular fingerprint regions in the MIR [ ] . Frequency 
combs [ , ] in the VIS to MIR enable ultrahigh resolu- 
tion spectroscopy used, for example, to study the compo- 
sition of extrasolar planets. Finally, when combined with 
advances in pulse shaping, coherent supercontinua enable 
arbitrary multi-color waveform synthesis of spectra span- 
ning the UV to MIR, which offers the possibility of simul- 
taneously probing and controlling [15] electronic, vibra- 
tional and rotational motion, or synthesizing single electric 
field waveforms at arbitrary carrier wavelengths [16]. 



These applications place widely varying demands on 
the characteristics of a supercontinuum (SC) source. 
Tremendous progress has been made using novel optical 
fibers, with SC spanning the UV into the MIR [5, 6, 17]. 
However, besides the intrinsic peak power limitations and 
alignment sensitivity, the most extreme bandwidths ob- 
served in fibers have been achieved at the expense of pulse 
to pulse repeatability, which often manifests in compli- 
cated temporal pulse splitting translating into spectral 
modulation [ , ]. A requirement for shot-to-shot re- 
peatability and a smooth spectrum therefore limits the 
allowable propagation length, restricting the bandwidth 
significantly [6]. 

A maturing understanding of the mechanisms of fil- 
amentation [19] reveals the difficulty of achieving signifi- 
cant spectral broadening to the long wavelength side of the 
pump: most of the underlying processes tend to blue shift 
the spectrum. One contribution to the spectral broad- 
ening is due to the free electrons which are produced by 
ionization throughout the pulse. These produce an ul- 
trafast decrease in the refractive index, leading to a blue 
shift. The effect may be clamped by ionization losses, 
dispersion, and defocusing. Another mechanism for spec- 
tral broadening is self-phase modulation, which produces 
a spectral shift proportional to the temporal intensity gra- 
dient, with red shift on the rising edge and blue on a 
falling edge. A key role, therefore, is played by tempo- 
ral reshaping of the pulse. Two of the main sources of 
temporal reshaping — self-steepening and space-time fo- 
cusing — intrinsically produce a steep trailing edge, again 
resulting in a blue shift. Other contributions to reshap- 
ing — ionization-induced defocusing and loss — also tend 
to act more strongly on the trailing edge and centre of 
the pulse respectively. Steep intensity gradients or shocks 
are naturally arrested by dispersion. Which of the mech- 
anisms dominate depends strongly on the pulse and ma- 
terial parameters [20, 21, 22]. Despite its pessimistic out- 
look for achieving significant spectral broadening to longer 
wavelengths, the modern elucidation of bulk material fil- 
amentation does raise the tantalizing prospect that, since 
the clamping of the blue shift is largely determined by 
the material properties, the short wavelength edge of the 
spectrum could be largely independent of the pump wave- 
length, offering the potential for enormous spectra when 
pumping in the MIR. More rigorous theoretical support 
for this conjecture comes from a Born approximation-like 
view of SC in bulk materials [ ] in which the spectral 
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Figure 1: (a) Measurement and (b) simulation of the far field {O^X} representation in absolute units; the grey area 
in (a) indicates the noise level and the detector ranges are shown above the figure, (c) Temporal intensity profile and 
(d) spectrogram of the on-axis simulated pulse emerging from the YAG plate; note the logarithmic intensity colour 
and wavelength scale in the spectrogram. 



extent is largely determined by the medium dispersion. 
In this view, the angularly-resolved ({6>, A}) spectrum is 
linked to the mechanism through phase-matching of "x- 
waves" . Normal, near-zero or anomalous dispersion man- 
ifests in "x", "fish", or "o" shaped structures [24, 25] in 
the {^, A} domain. 

Here we present the broadest supercontinuum ever pro- 
duced in bulk, from filamentation of a femtosecond MIR 
pulse in a thin plate of yttrium aluminium garnet (YAG). 
We used our homebuilt optical parametric chirped pulse 
amplifier (OPCPA), which delivered CEP stable, 10 /iJ en- 
ergy, 85 fs pulses at 3100 nm at 160 kHz repetition rate [ ]. 
6.9 fiJ of the output was collimated to 6.8 mm and focused 
by a 75 mm focal length CaF2 lens to a spot size of 50 /im. 
The peak power was 76 MW which is three times above 
the critical power for YAG. Beam sizes are given as ra- 
dius at the 1/e^ value of the peak intensity. The 2 mm 
thick YAG plate was placed into the focal plane (i.e. the 
front surface is located 1mm from the focal plane), pro- 
ducing a filament inside the YAG plate which is clearly 
visible to the naked eye. We measured a {^, A} spectrum, 
the main challenge being the tremendously large wave- 
length range. We covered the spectrum by simultaneously 
using three intensity-calibrated fiber-coupled spectrome- 
ters. They consisted of a silicon CCD spectrometer for 
wavelengths up to 925 nm, an InGaAs CCD spectrometer 
for 925-2500 nm, and a Fourier-transform infrared spec- 
trometer (FTIR) with MCdTe detector for wavelengths 
above 2500 nm. The three fibers were held parallel to each 
other with their tips separated horizontally by roughly a 
millimeter. They were then scanned vertically across the 
beam, with spectra being acquired at each position. An 
image registration step in the processing accounted for the 
different positions of the fibre tips. 

Figure 1(a) shows a measured {6>, A} spectrum. The 



long wavelength edge decays smoothly, reaching the noise 
level at 4500 nm, whilst the short wavelength edge has a 
sharp peak at 610 nm before dropping rapidly into the 
noise at 400 nm. Several structures are present which 
we analyze below. To obtain higher sensitivity and also 
demonstrate that all spectral components of the SC may 
be practically collected and relayed, we refocused the SC 
into each spectrometer fibre with a 50 mm lens. The re- 
sulting spatially averaged spectrum is shown in Fig. 2. 
This measurement featured a lower pulse energy but a 
tighter focus, leading to a slightly different spectrum, and 
we used different detectors above 2500 nm as indicated. 



The highest spectral energy density is found around 
the MIR pump wavelength with a maximum energy con- 
tent of lOnJ/nm. The smallest energy content is ob- 
served in the 750-1000 nm range at a few pJ/nm; note 
that these energies are still sufficient for most absorption 
spectroscopy experiments. Figure 1(a) clearly shows a 
"fish-tail" structure at 610 nm in the {^, A} representation 
whose energy content is roughly one order of magnitude 
higher than at the minimum at 850 nm. A comparison 
with theory shows that this feature is not related to odd 
harmonics of the MIR driver; rather it is consistent with 
the aforementioned dynamics of x- waves. A small peak 
can also be observed at around 1100 nm, which simula- 
tions show contains a strong contribution from the third 
harmonic of the pump. 
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Figure 2: Supercontinuum generated by 3100 nm, 2.6 /iJ 
pulses in YAG. The ranges of each spectrometer /detector 
are indicated, as weh as the noise level for the long wave- 
length detector (grey horizontal line). 

A practical issue for applications of a SC is the shot- 
to-shot reproducibility and coherence across the spectrum. 
The process must therefore be deterministic and not exces- 
sively sensitive to the inevitable shot-to-shot fluctuations 
of the drive laser. These properties have been demon- 
strated in bulk supercontinuum generation with 800 nm 
pumping [27, 28, 29], and we now show that they also 
hold for our MIR pumping conditions using an f-2f in- 
terferometer. In the interferometer, the filament was re- 
focused into a type I beta barium borate (BBO) crystal 
resulting in second harmonic generation from a wavelength 
range of ^ (1500 ± 300) nm. The frequency-doubled com- 
ponent interfered with the short-wavelength edge of the 
SC at ~ (750 ± 150) nm, producing spectral fringes ac- 
quired with a spectrometer. Figure 3 shows a 15 s subset 
of such a dataset recorded over a 10 minute interval. The 
amplitude of the observed fluctuations provides an upper 
bound on both of the potential sources of phase jitter. 
Over the full 10 minutes we obtain a standard deviation 
of 259mrad. 
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Figure 3: Spectral interference fringes produced by super- 
continuum in / — 2/ interferometer versus time. 

Figure 1(b) shows the A} spectrum obtained from 
numerical simulations of the nonlinear propagation of the 
IR laser pulse in the 2 mm long YAG plate. A comparison 
with Fig. 1(a) shows very good qualitative agreement be- 
tween measurements and simulations. The most intense 
part of the spectrum corresponds to the "o"-like structure 
typical of the anomalous dispersion regime. The hyper- 
bolic branch crossing the axis around 560-620 nm corre- 
sponds to the tail of the "fish" -wave, i.e., the loci of phase 
matched angle and wavelength lying in the visible region 
of the spectrum. Simulations show that this tail appears 



in the last 400 /im of the YAG plate and is not affected by 
third harmonic generation (THG). Switching off THG in 
simulations by e.g. propagating the envelope rather than 
the carrier resolved field only suppress the peak emerging 
around 1000 nm. This also indicates that CEP does not 
affect the "fish" wave structure of the spectrum. These 
features are always present and start with a slightly dif- 
ferent focal position with respect to the entrance face of 
the YAG plate. They also appear with different Kerr coef- 
ficients and ionization rates, the accuracy of which is not 
known with great certainty. Increasing n2 up to a factor 
of 4 leads to a broader extent of the main spectral peak 
centered at 3000 nm and a change in the balance between 
THG and the SC. Higher Kerr coefficients lead to a pre- 
vailing "fish" tail with respect to the third harmonic. 

The good agreement between experiment and simu- 
lation permits a deeper study of the temporal dynamics 
using the latter. Figure 1(c) and (d) show the emerging 
on-axis temporal pulse profiles. Figure 1(c) shows that 
pulse splitting has not occurred over the nonlinear propa- 
gation distance since the main pulse has developed a pre- 
pulse which is not split off completely from the main pulse. 
These two parts are separated by only a few femtoseconds 
and preceded by a shallow increasing background. A more 
detailed picture is revealed in the spectrogram (Fig. 1(d)) 
with a logarithmic intensity color scaling. We find that the 
short wavelength extension — ending in the "fish tail" — 
consists of a low-intensity trailing edge at 25 fs. The very 
weak spectral modulation observed in the simulated {^, A} 
representation between 1000-2000 nm is attributed to the 
interference of the main pulse with its trailing pedestal. 
The amplitude of this pedestal is however critically de- 
pendent on propagation length and uncertain nonlinear 
parameters such as the nonlinear refractive index and the 
ionization rate. We clearly do not observe these weak mod- 
ulations in our measurement; we attribute this to an ab- 
sent or much less pronounced pre-pulse or pedestal. 

In this Letter we demonstrated the first stable multi- 
octave SC from filament ation of MIR femtosecond pulses 
in bulk material. We measured a smooth spectrum from 
450 nm to 4500 nm, corresponding to 3.3 octaves, with 
a spectral energy density of 2pJ/nm to lOnJ/nm. The 
smoothness of the spectrum indicates absence of chaotic 
and multiple pulse splitting and a CEP measurement indi- 
cated good shot-to-shot coherence over a significant spec- 
tral range. We obtained an angularly-resolved far-field 
spectrum which shows good agreement with full 3D simu- 
lations. These indicate that the mechanism is dominated 
for formation of an "o" wave, with the short wavelength 
edge described by phase-matching. By demonstrating a 
simple and robust method for coherently extending the 
spectrum of an amplified femtosecond MIR pulse to the 
UV we provide a powerful platform for the next genera- 
tion of ultrafast applications. 
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